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Abstract A three-dimensional finite solid element model for 
Japanese traditional post-and-beam connections was devel¬ 
oped using the wood foundation method, which employed the 
concept of a beam on a nonlinear foundation. The wood foun¬ 
dation in the model was a three-dimensionally prescribed zone 
surrounding a nail shank in order to address the intricate wood 
crushing behavior induced by nail slip. Material models for the 
wood members and the foundation were developed based on 
the transversely isotropic plasticity from the software package 
ANSYS. The Japanese post-and-beam connection modeled 
was a ten-nail multiple connection with a mortise and tenon 
joint and is called the CPT (Comer Plate, T-shaped) connec¬ 
tion. Details of the model development are presented. As a 
feasibility study, blind predictions of the model were compared 
with available connection test data and showed good results 
for predicting the progress of the load-deformation relations 
in three dimensions. However, the limitation of the model was 
found in simulating fracture failures such as wood splitting or 
nail tear-out from the wood. Model applications and the need 
for model improvement are discussed. 

Key words Three-dimensional finite element model ■ Nailed 
connection ■ Wood mechanical properties 


Introduction 

Utilization of Canadian wood species for Japanese tradi¬ 
tional post-and-beam construction has been comprehen- 
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sively studied over past 10 years at the University of British 
Columbia (UBC), Canada. 12 As a part of the studies, the 
development of three-dimensional finite element (FE) 
models for the Japanese traditional post-and-beam connec¬ 
tions was initiated. The objective of the modeling work was 
to develop deeper understanding on three-dimensional 
(3D) stress states in a complicated connection system in an 
effort to develop innovative solutions for improving the 
Japanese traditional connection method with Canadian 
wood products. Since the traditional post-and-beam con¬ 
nection includes the interaction of the mortise and tenon 
joint and the nailed connections, FE analysis with solid ele¬ 
ments was considered the best approach to investigate the 
behavior of the connection under loads. 

However, the lack of current techniques in finite solid 
element modeling for wood materials in connections neces¬ 
sitated an original study to develop a 3D FE model for 
nailed connections because the Japanese post-and-beam 
connection includes multiple-nail connections with a 
T-shaped side metal plate. 

To date, 3D FE models of wood connections have focused 
mainly on bolted connections because the crushing (or com¬ 
pression) behavior of wood under a bolt with a relatively 
large diameter was assumed to follow the conventional 
mechanical properties of wood that were determined 
through uniaxial compression tests. 3-7 However, this assump¬ 
tion is not applicable for nailed connections because the 
interaction between the nail and the wood is governed pri¬ 
marily by the crushing behavior of the local wood contact¬ 
ing the nail rather than by the gross compression behavior 
of the wood. 1. 

Recently, a 3D FE modeling method that is capable of 
accounting for wood crushing behavior, called the 3D wood 
foundation model, was introduced by Hong et al. 910 A 
model using a 3D wood foundation for single-nailed con¬ 
nections showed good agreement in the ramp load-defor¬ 
mation curves and the deformed shapes. In this study, 
therefore, a feasibility study of the finite solid element 
model for multiple-nailed connections was conducted with 
the Japanese traditional post-and-beam connection, called 
the CPT (Comer Plate, T-shaped) connection. Blind predic- 
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tions of the model were compared with the test data of the 
CPT connection that were available from a previous test 
program at the University of British Columbia. 1,11 All rele¬ 
vant FE analyses in this study were performed using the 
commercial software package, ANSYS 


parameters were assigned to the prescribed foundation 
zone in order to account for the load-wood crushing 
(embedment) behavior. 


Material parameters for the 3D FE models for the wood 
and wood foundation 


Theory 

Three-dimensional wood foundation model in dowel-type 
wood connections 

Foschi and Bonac 8 developed the wood foundation method 
for modeling a dowel-type wood connection that treats the 
wood under the connector as a nonlinear foundation. In this 
study, the three-dimensional wood foundation method 
introduced by Hong et al. 9,10 was used for modeling the 
nailed connections of the CPT connection. The three- 
dimensional wood foundation is a prescribed zone sur¬ 
rounding a nail shank. Figure 1 shows a half-size analog of 
the wood foundation used for a nail-embedment solid 
element model. The foundation material parameters were 
determined through nail-embedment tests, and then the 



Fig. 1 . Conceptual illustration for the wood foundation model in the 
three-dimensional finite element nail-embedment model, d is the diam¬ 
eter of the nail, M is the multiplier for determining the depth (M x d) 
of the wood foundation 


The procedures for determining the material parameters 
used for the 3D wood foundation method are introduced in 
Table 1. The transversely isotropic plasticity model with 
bilinear stress-strain relations, as used by Moses, 4 ' 5 was 
adopted for 3D FE material modeling of the wood and 
wood foundation. This model was implemented through the 
TB,ANlSO option in ANSYS, which is based on Hill's 
potential theory and an associative flow rule. 4512 Based on 
compressive strength properties, equality of bilinear stress- 
strain relationships in tension and compression was 
assumed. 4-7 Transversely averaged Poisson’s ratios taken 
from the literature were used. 13 Yield points in the bilinear 
stress-strain relationships were determined by bilinear 
curve fitting to the load-deformation curves from wood 
compression tests. 

For the wood foundation, the bilinear stress-strain rela¬ 
tionships were derived from the nail-embedment curves 
parallel and perpendicular to the grain direction as illus¬ 
trated in Fig. 2. The nominal bilinear stress-strain relations 
for the foundation were calculated by: 910 





A 

w y 


( 3 ) 


where a nom is the nominal yield stress of the wood founda¬ 
tion (MPa); e nom is the nominal yield strain of wood founda¬ 
tion (mm/mm); K is the nominal foundation modulus (MPa); 
d is the diameter of the dowel (nail) shank (mm); P y is the 
yield load in the bilinear load/unit length-embedment curve 
(N/mm); and W y is the yield deformation in the load/unit 


Table 1. Procedures for determining the material parameters of the three-dimensional (3D) finite element (FE) models for the wood and wood 
foundation 


Requisite constant 

Symbol 

Direction 

Method of determination 

Initial modulus 

E 

L, R, T 

Compression test, or dowel embedment test 

Initial shear modulus 

G 

RL, LT, RT 

Coupled bilinear constitutive model a 

Poisson’s ratio 

V 

RL, LT, RT 

Average values from the literature b 

Tensile yield strain 

e y 

L, R, T 

Equal to compressive properties 

Tensile tangent modulus 

E 

L, R, T 

Equal to compressive properties 

Compressive yield strain 


L, R, T 

Compression test, or dowel embedment test 

Compressive tangent modulus 

E 

L, R, T 

0.01 x initial modulus ( E ) 

Shear yield strain 

7, 

RL, LT, RT 

Coupled bilinear constitutive model a 

Shear tangent modulus 

G' 

RL, LT, RT 

0.01 x initial shear modulus (G) 


The parameters in R and T directions are equal according to the transversely isotropic plasticity 

L, longitudinal; R, radial; T, tangential 

“The theoretical models developed by Saliklis et al. 15 

b In this study, the values from the Wood Handbook 13 were cited, except for the computed Poisson’s ratio as v TL = v LX • (E X /£ L ) 
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Embedment (mm) 


Fig. 2. Example of a load per unit length-embedment curve 


length-embedment curve (mm). For the nominal yield 
stress and strain, the diameter of the nail shank (d) was 
chosen to normalize the dimensions of variables involved 
since the diameter of the nail shank was presumed to be the 
most important nonmaterial factor governing the localized 
wood crushing behavior. Note that Eq. 1 for the nominal 
yield stress is the same form as given for calculation of 
dowel-embedment strength in ASTM D 5764. 14 Once the 
nominal foundation properties are determined, these should 
be calibrated to become the effective properties depending 
on the geometry of the foundation. For shear strength prop¬ 
erties of the wood and wood foundation, shear moduli and 
shear yield points were estimated using the theoretical rela¬ 
tionship given by Saliklis et al.: 15 


Initial shear modulus: G, ( = G jt = 


Shear yield strain: (y y ) = (y y ) = 


jEfE, 

(4) 

2-(l + Vvv„.) 

(<7int.)i 

(5) 

1.98-Ej \ G tj 


where i, j are the parallel-to-grain (L) or perpendicular-to- 
grain (T) direction; E is the initial modulus (compressive 
modulus of elasticity, or effective foundation modulus 
(MPa); G is the initial shear modulus) (MPa); v is Poisson’s 
ratio; y y is the shear yield strain in the bilinear shear stress- 
strain curve (mm/mm); and cr int is the stress intercept in the 
bilinear normal stress-strain curve (MPa). The tangent 
moduli in compression and shear were taken empirically as 
a value of 0.01 x initial modulus. Depending on the dimen¬ 
sion of the prescribed foundation, the nominal foundation 
modulus ( K ) and yield strain (e nom ) should be calibrated to 
the effective foundation properties in initial modulus ( E ) 
and yield strain (el), as shown in Table l.The effective foun¬ 
dation properties were determined according the calibra¬ 
tion procedures given by Hong et al. 910 These effective 
foundation properties should enable the 3D FE nail-embed¬ 
ment model to simulate the load-embedment curves in par¬ 
allel- and perpendicular-to-grain directions simultaneously 
(e.g., see Fig. 3). 



Fig. 3. Simulated and tested load-nail embedment curves in parallel- 
to-grain direction for a nail diameter of 3.3 mm and a wood specimen 
38 mm thick. The no-foundation model used only conventional wood 
material parameters derived from uniaxial compression tests. This 
resulted in an unacceptably stiff curve 
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Fig. 4. An example of the Japanese post-and-beam connection test 
carried out in a University of British Columbia research program 
(units: mm) 


Materials and methods 

The Japanese CPT connection is a standard connection for 
attaching a post to a beam in Japanese traditional post-and- 
beam construction. This connection consists of a T-shaped 
steel plate connector, called a CPT metal plate, with ten 
Japanese standard nails. Often, a mortise and tenon joint is 
combined to strengthen the connection. Figure 4 shows the 
CPT connection and the test setup used for the tension test 
in the UBC research program. 111 The test configuration and 
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Table 2. Material parameters of 3D FE model using Douglas fir 



Initial modulus, 

Initial modulus, 

Yield stress 

Yield stress 

Poisson’s ratio 

Poisson’s ratio 



E r = E t 

o L 

°r = O t 

V LT 

v RT 


(MPa) 

(MPa) 

(MPa) 

(MPa) 

() 

o 

Wood 

16900 

832 

43.9 

5.03 

0.37 

0.38 

Wood foundation 

736 

142 

16.0 

6.53 

0.37 

0.38 


L, R, and T represent longitudinal, radial and tangential directions, respectively 


Fig. 5. Three-dimensional finite 
element CPT (Comer Plate, 
T-shaped) connection model (a) 
and the boundary conditions (b) 



the results were used for the blind prediction of the 3D FE 
model. A CPT metal plate was installed with ten Japanese 
standard ZN65 nails (hand driven) on one side of a post- 
and-beam member with a 105 mm x 105 mm cross section. 
The wood members of the UBC connection specimens were 
Canadian coastal Western Hemlock (Tsugci heterophylla), 
which had a range of specific gravities (SG) from 0.39 to 
0.47. In this feasibility study, however, the CPT connection 
model was developed with Douglas fir (Pseudotsuga men- 
ziesii , average SG = 0.54), on which the development of the 
original wood foundation model had been based. 9 ' 10 Further 
studies should be considered for the generalized use of the 
model with other wood species. Douglas fir wood, con¬ 
ditioned to a moisture content of 15%, was used for the 
material property tests. According to ASTM D 143, the 
compression tests were conducted using a 25-mm wide x 
25-mm thick x 100-mm high wood specimen in the longitu¬ 
dinal, radial, and tangential direction, respectively, with 45 
repetitions. Also, for the material properties of the wood 
foundation, nail-embedment tests were conducted using 
Japanese standard nail ZN65, which has a shank diameter 
of 3.3-mm and a length of 65-mm. The tests were done in 
longitudinal, radial, and tangential directions with 15 repeti¬ 
tions, respectively, according to the ASTM D 5764 half-hole 
embedment test. The final material parameters for the 
model input of the wood and wood foundation are sum¬ 
marized in Tabic 2. For the material model of nails, an 
elasto-perfectly plastic model was assumed with a yield 
stress of 360 MPa, a modulus of 200 GPa, and a Poisson’s 
ratio of O.3. 9 For the CPT plate connector, a yield stress of 
250 MPa, a Young’s modulus of 200 GPa, and a Poisson’s 
ratio of 0.3 were used. 111 


Finite solid element modeling 

The dimensions and boundary conditions of the CPT connec¬ 
tion model followed the UBC test configuration as shown in 
Fig. 4. No influence from loose or tight wood-to-wood contact 
in the mortise and tenon joint was considered. 

To embody the connection, eight-noded, quadrilateral iso¬ 
parametric brick elements (SOLID45 from ANSYS) were 
used. The surface-to-surface contact elements were defined 
on every contact interface with a frictional coefficient of 0.7 
for wood-to-wood and wood-to-steel contacts, and 0.3 for 
steel-to-steel contact (CONTA174 and TARGE170 from 
ANSYS). 4,16 Figure 5 shows the solid element model with the 
post-and-beam member, the nails, and the plate connector 
separated, and a description of boundary conditions. 

The size of the wood foundation was chosen to have a 
square cross section with a depth of 4.5 x 3.3-mm nail diam¬ 
eter (M x d as shown in Fig. lj.The selection of the multi¬ 
plier as 4.5 was based on the efficiency of solution 
convergence for the simulation of nail-embedment behav¬ 
ior studied by Hong et al. 910 Since the nail spacing of the 
connection was narrow, two adjacent foundations over¬ 
lapped partly. A wood foundation with a square cross 
section facilitated the meshing on the overlapping regions. 


Results and discussion 

Deformation and failures of the CPT connection 

Figure 6 shows the progress of the simulated deformation 
when displacement-controlled, and force (pressure )-con- 
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Fig. 6. Simulated displacement 
(sectioned views) of the CPT 
connection on a displacement- 
controlled loading and b 
force-controlled loading 



2.0 mm 3.6 mm 7.0 mm 

(a) displacement-controlled loading 





5.8 kN 



(b) force-controlled loading 



Fig. 7. Actual failure modes of 
the CPT connection under 
tension: a splitting on tension 
perpendicular to the grain, b nail 
pull-out and CPT plate in shear 
on a specimen with a short 
tenon 1,2 



(a) 


(b) 


trolled loading were applied. The displacement-controlled 
loading has all nodes on the top surface of the post uni¬ 
formly displaced in the direction parallel to the post, while 
in the force-controlled loading, the applied forces always 
acts in a direction normal to the top surface of the post, 
wherever the post moves. 

The simulated deformation using displacement-con¬ 
trolled loading showed better agreement with the actual 
deformation shown in Fig. 7. The difference resulting from 
the different loading methods in the simulated deformation 
was the tilting behavior of the post member. The model 
using incremental displacement-controlled loading showed 
a uniformly vertical translation of the post member while 
the force (pressure)-controlled loading represented the 
tilting behavior (see Fig. 6b). This tilting behavior did not 
agree with the actual deformation. The simulated results 
using the force-controlled loading were, therefore, aban¬ 
doned for analysis purposes. 


Figure 8 shows the simulated deformation of nails and 
CPT plate under displacement-controlled loading. The sim¬ 
ulated flexural deformation of the three nails penetrating 
the tenon agreed well with the observed nail pull-out from 
the bottom sill and intensively plastic bending shown in Fig. 
7b. Also, the simulated plastic strain contours in Fig. 8b show 
a good indication of the progress of the plate shear failure. 

Modified model 

The major failure modes of the CPT connection were per- 
pendicular-to-grain tension splitting of the beam member, 
the nail pull-out, and the plate shear. However, during the 
process of model development, it was discovered that the 
most influential failure on model prediction was the nail 
tear-out from the tenon as, shown in Fig. 9. 

In the CPT connection, the three middle nails of the five 
nails in the beam member penetrated the tenon. In this 
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Fig. 8. The Y-directional plastic 
strain contours in the nails (a), 
and the CPT plate (b) 
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Fig. 9. Typical failure modes of the three nails penetrating the tenon: 
nail tear-out from tenon 



Fig. 10. Simulated deformed shape of the modified CPT connection 
model 


configuration, the three nails had a narrow end or edge 
distance in the tenon. Therefore, it is most likely that the 
failures of the nail connections within the tenon occurred 
in a brittle manner, including the edge splitting of the side 
nails and the tear-out of the centre nail. These kinds of 
failures may occur at the initial deformation level, or even 
when a nail is driven improperly. Once a nail is pulled out 
of the tenon, the load-carrying ability of the nail connection 
within the tenon vanished. The model, however, could not 
simulate this kind of fracture process. This is a limitation of 
the model in this study. As shown in Fig. 10, the wood ele¬ 
ments surrounding the centre nail could not split off unless 
a fracture algorithm or an intentional cleavage for allowing 
the tear-out was introduced into the model. As long as the 
nail was surrounded by wood elements, the load-carrying 
ability of the connection was retained in the model. In order 
to investigate the influence of these initial fractures on the 
simulation, the original CPT connection model was modi¬ 
fied to remove the load-carrying ability of the side nail 
connections in the tenon, by eliminating the contact ele¬ 


ments between the side nails and the tenon. The nails could 
then move independently, regardless of the surrounding 
wood elements in the tenon. The centre nail connection 
model remained without any modifications. 

Simulated load-deformation curves 

The load-deformation curves simulated using the original 
and modified connection models are superimposed on the 
test data in Fig. 11. Note that the two black bold curves in 
Fig. 11 are the test data of the CPT connection with a short 
tenon that the three middle nails did not penetrate (see the 
tenon in Fig. 7b). In the tests and the simulation, the three 
middle nail penetrations seem to significantly influence the 
initial stiffness or the load-carrying capacity of the connec¬ 
tion, although the existence of the initial fracture in the real 
tenon was unknown. This implies that if 3D FE parametric 
and stress-quantitative studies are conducted for optimizing 
the configuration of the CPT connection, some improve¬ 
ment for the connection performance should be explored. 
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Fig. 11. Simulated load-deformation curves of the initial and the modi¬ 
fied CPT connection model 


Conclusions 

A three-dimensional finite element model of a Japanese 
traditional post-and-beam connection with multiple-nail 
connections and a wood-to-wood joint was developed using 
the three-dimensional wood foundation approach. This 
three-dimensional model showed a good ability to reveal 
the internal behavior of the connection. Simulated results 
that can be quantified three-dimensionally in terms of stress 
and strain would lend themselves to the improvement of 
dowel-type wood connections. 
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